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The purpose of this study was to investigate the interaction of 
plane elastic waves with a thin, hollow, cylindrical shell embedded in an 
elastic medium. 

Thetevsimemicaleonell 1s considered) to be elastic, isotropic, 
homogeneous, and of infinite length. It is surrounded by an elastic, iso- 
tropic, and homogeneous medium whose motions conform to the ordinary theory 
of elasticity. A plane stress wave, either dilatational or shear with a 
step variation in time, whose wave front travels in a direction perpendic- 
ular to the cylinder axis, envelops the shell. Later, a Duhamel integral 
is used to study other wave shapes for the incident stress. 

The response of the shell is studied by expressing the two com- 
ponents of displacement, radial and tangential, in terms of Fourier series, 
each term of which is called a mode. The equations of motion of the shell 
in vacuo are derived from expressions giving the strain and kinetic ener- 
gies due to generalized external forces. Forces on the shell result from 
the stresses in the medium at the boundary. Stresses in the medium are 
taken to be the sum of the stresses due to the incoming stress wave expressed 
im Germs Of Bourmiersserics whose coefficients are known, and those due to 
the reflected and diffracted effects expressed in terms of a pair of dis- 
placement potentials representing waves diverging from the axis of the 
shell, 

Inewequattons to be solved consist of two pairs of coupled integro- 


differential equations in the generalized coordinates of the shell and the 





displacement potentials. By use of a digital computer they are solved mode- 
wise by a step-by-step iterative integration technique known as the Newmark 
Beta Method, with which values for the potential functions, and the accelera- 
tions, velocities, and displacements of the shell are determined. Stresses 
in the shell are found from the displacements, and the values of the potential 
functions permit determination of stresses for any point in the medium. 

Although the equations are written to include an infinite number 
of modes, only the first three modes are considered in detail. The computed 
solution is compared with values obtained from a series expansion of the 
equations, which is valid for short times, and with the static solution based 
on the theory of elasticity to which the general solution should approach 
asymptotically. In addition, the results of two particular problems are 
compared with results given in another study. 

Numerical solutions are obtained to determine the effect of the 
several parameters which describe the relative physical properties of the 


shell and medium. Results are presented in tabular and graphical form. 
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CHAPTER 1 


INTRODUCTION 


1.1 General Remarks 

The problem of designing underground protective structures to 
resist the effects of nuclear weapons has become increasingly important in 
recent years with the development of modern weapons whose destructive capacity 
is overwhelming. Engineers in this field are hampered to a great extent by a 
lack of theoretical information on how structures in media such as soil or 
rock behave when subjected to dynamic loads. Even for static loads alone, 
much of the design practice today is of a semi-empirical nature. 

When a nuclear explosion occurs, stress waves are transmitted 
through the air and ground. How are they transmitted and how are they modified 
by the presence of a structure embedded in the medium? How does the structure 
respond? 

The purpose of this report is to study one aspect of the problem, 


the interaction between the medium and structure. 


1.2 Statement of Problem 

The problem considered here consists of analyzing the elastic 
response of a hollow cylindrical shell (tunnel lining) embedded in an elastic 
medium when subjected to an incident plane stress wave traveling in a direction 
perpendicular to the axis of the shell. Some questions with which this problem 
may be associated are: Do tunnel linings in contact with rock afford a measure 
of protection significantly higher than an unlined opening? What magnitude 
AidmbiMmesvariavlonrOlsdisplacenenu,. velocity, and acceleration would equipment 


mounted within such a structure be subjected to? How are stress waves within 





the medium modified in the vicinity of the shell? This study was conducted 
in an attempt to find some qualitative and quantitative answers to these 


questions within the limitations imposed by the assumptions noted below. 


1.3 Basic Assumptions 


The cylindrical shell is considered to be of infinite length, and 
is embedded in an elastic medium of infinite extent in all directions. A 
plane stress wave whose front travels in a direction perpendicular to the 
cylinder axis envelops the shell. Strains parallel to the axis in both the 
medium and shell are assumed to vanish; thus, since each cross section of the 
shell is exactly similar to every other, the problem is reduced to one of 
plane strain. 

In the mathematical development of the problem certain basic 
assumptions were made, the most important of which are given here with a few 
explanatory remarks: 

(1) The medium is considered to be homogeneous, isotropic, and 
linearly elastic. This implies that the ordinary theory of stress wave 
propagation applies. In view of the non-homogeneous, non-isotropic, and 
non-elastic characteristics of most materials encountered in nature, this 
is a severe limitation; however, current theories of stress propagation 
through such media have not advanced to the stage where this limitation can 
be readily overcome. In the case of some rocks, though, this assumption 
may be justified. 

(2) The material in the shell is also considered to be homogeneous, 
isotropic, and linearly elastic. Generally speaking, this assumption is valid 
for values of stress below the so-called proportional limit of materials 


Commonly sused.. in addition, the thickness of the shell relative to its radius 





is assumed small; this permits expression of all stress components of the 
shell in terms of a function which describes the deflection of its middle 
surface. This deflection must satisfy a linear partial differential equation 
with the appropriate boundary conditions. 

(3) The incident stress wave considered is either a plane dilata- 
tional or a plane distortional (shear) wave. Under actual conditions, both 
waves are propagated with the shear lagging the dilatational wave. The 
combined effect for elastic conditions may be determined through the principle 
of superposition. 

(4) The radial and tangential particle velocities of the medium 
at the boundary are equal to that of the shell. This is the continuity 
relation insuring that the shell and the medium are in contact with no 
relative slip occurring at the boundary. 

(5) Any additional mass within the shell is assumed to be dis- 
tributed symmetrically about the axis. The significance of this assumption 
is found in the development of the equations of motion to account for any 
additional mass located within the shell. 

Other assumptions are presented in the formal development of the 


mathematical expressions used to describe the behavior of shell and medium. 


1.4 Method of Approach 

The two components of shell displacement, radial and tangential, 
are written in terms of Fourier series from which expressions giving the 
strain energy and kinetic energy of the shell in vacuo are derived. The 
equations of motion are derived from Lagrange's equations in terms of the 
displacement functions and forces acting on the shell. Forces on the shell 


result from the stresses in the medium at the boundary. Stresses in the 





medium are taken to be the sum of the stresses due to the incoming stress wave 
expressed in terms of Fourier series whose coefficients are known, and those 
due to the reflected and diffracted effects expressed in terms of a pair of 
displacement potentials representing waves diverging from the axis of the 
shell. The form of these potentials as derived by Lamb Gee 1s in terms of 
sine and cosine series, the nature and treatment of which has been studied 
iby Paul (6). 

The equations of motion, described as a pair of coupled integro- 
differential equations, are solved modewise using a numerical technique 
known as the Newmark Beta Method (5) which permits determination of the 
coefficients of the potentials and values of acceleration, velocity, and 
displacement of the shell. 

The solution obtained is compared with values obtained from a series 
expansion of the equations, which is valid for short times, and with the static 
solution based on the theory of elasticity which the machine solution should 


approach asymptotically. 


1.5 Previous Work 

The problem stated above has been the subject of a recent report 
by Baron (1), whose analysis consists of first solving for the displacements 
caused at the boundary of an unlined cylindrical cavity subjected to a plane 
stress wave. This is done through an integral transform approach, the solu- 
tion of the transformed equations being expressed using Hankel functions. 
The evaluation of the inverse transform is accomplished only with great 
computational etieort. | Values of displacements obtained are then used as 


influence coefficients in determining the displacements of the shell. 





* Numbers in parenthesis refer to the corresponding entry in the Bibliography. 





Solutions cbtained for two different sets of parameters by Baron are compared 
in Section 4.6with results obtained by the method of solution outlined in 
this Teporu. 

The study by Paul (6) consisted of analyzing the effect of a plane 
stress wave incident on an unlined cylindrical cavity in an elastic medium. 
The reflected and diffracted waves are described in terms of displacement 
potentials which represent outgoing shear and dilatational waves. <A method 
was developed for determining values of these potentials, and a similar 


method is used in this report. 


1.6 Notation 

Notation is defined throughout the text where it first appears; 
however, the following list summarizes the main uses of certain symbols. In 
discussions of special topics other meanings may be ascribed to the symbols, 


at which time they will be redefined. 


A = Area of cross section of the shell, per unit length. 

A,» B C. = Coefficients of Fourier series for stresses in the 
medium. 

AL B Cc = Coefficients of Fourier series for total stresses 


in the medium. 


AM 2 EM =  Weiehting factors. 

a3 ice = Generalized coordinates for the displacements of 
the medium at the boundary. 

a,b = Generalized coordinates for displacements of the 

ns ns 

shell. 

Cy» Cy = Velocities of wave propagation, dilatational and 
shear, respectively. 

d = Distance trem Neutral axis of shell to its outermost 


iupers 





U 


& | 


Moduli of elasticity for the medium and shell, 
respectively. Bar over the symbol refers to the 
plane strain modulus. 


Volumetri eG stra in. 


Generalized coordinates of the displacement 
potentials. 


MOnenGiol inertia of che shell per init. veneain 


Parameter which relates to the shape of a time 
dependent stress wave. 


By gee ac ote mea cael 

alone ratio of velocities. 

Mass of the shell per unit surface area. 
Additional mass within the shell. 

Mode number. 

Generalized forces. 

Radius of the shell. 

Maleiplyaines Tacters. 

Polar coordinates. 

Kinetic energy of the shell. 

Kinetic energy of additional mass within the shell. 
Thickness of the sheil; also, time. 
otrain energy of the shell. 

Displacement vector. 


Components of the displacement vector in the radial 
and tangential directions, respectively. 


Radial and tangential displacement components of 
the shell. 


Components of the displacement vector in the x and y 
directions, respectively. 


Rectangular coordinates. 





non-dimensionalized radial component of 
displacement. 


non=-dimensionalized tangential component 
p of displacement. 


Strains along the x and y axes, respectively. 
Circumferential strain. 

A variable of integration. 

Me eo bameue re 


PeoeboraveLer. 


peo pobanieter . 
» & parameter. 


7 a. Dabanletver. 


ele ro} ct oe I ‘| > ol ae 


(1+v) (1 -2y) 
dem) 


» & parameter. 

Polar cCeoerdinaves. 

Position angle of the incident wave. 
Curvature of the shell. 


Lame constants. 


Poisson's ratic. 


Vv 


Ley 
Arguments of the functions F and G, respectively. 
Mass density of the medium. 


Mass density of the shell. 


Normal components of stress parallel to the 
x and y axis, respectively. 


ohearing stress component. 





Lr 


01.9 
sb 


0 
smax 


g¢ 


Radial and tangential normal stresses in polar 
coordinates. Additional subscripts s and m, 
when used, refer“vo the shell and medium. 
Shearing stress in polar coordinates. 

Beading (stress 2a) ene snelie 

Maximum stress in the shell. 

Amplitude of incident dilatational wave. 
Amplitude of incident shear wave. 


Non-dimensionalized time in the case of the 
incident dilatational wave. 


Non-dimensionalized time in the case of the 
incident shear wave. 


Displacement potentials, dilatational and shear, 
respectively. 





CHAPTER IT 


BASIC EQUATIONS 


2.1 Equations for the Medium 


The differential equation of motion of a particle in a homogeneous, 
isotropic, and linearly-elastic medium in terms of its displacement vector u 
is given by Kolsky (3) in the form 


2 


Q 
c | 


CoP 2Unyiy = - uv x ve a) =— (2-1) 


ND 


ot 
where ’ and uw are the Lame constants, and p is the density of the medium. The 
vector u may be expressed as the sum of two displacements, the aveadient of 


a scalar potential and the curl of a vector potential 
u=Vo+t x ¥ (2-2) 


Here ) is a potential giving rise to an irrotational displacement and Ya 
potential leading to an equivoluminal displacement. 
If the wave equations (2-3) are satisfied and u expressed as in 


(2-2), the equations of motion (2-1) are automatically satisfied. 


a 
aC 
Po-528 
c, ot 
(22) 
2 1 OF 
Vi ete es 
a She 
a 
where the velocities of wave propagation are 
a wee 
1 p 
C2 
a ee 
C,, Sie 


for the dilatational and shear waves ,-respectively . 





to 


Since the problem is essentially a two-dimensional one wherein only the 
rotation about the cylinder axis is considered, the expression for the 
propagation of shear waves can be written in terms of a scalar potential 
function 


za — 
2 Ba Way Bet (2-3a) 
t 


<J 

ow 

I 
rs pole 


pm and y are functions of x, y, and t. The components of the displacement 


vector u can be expressed as 


n= 2+ (2-5) 


Stress components are 


aoe =e) CUE, 

GG  =py 2=6 
pea Te (2-6) 
O ee ee 

va J. 


where e represents volumetric strain, cy and Sa strain along the x and y axis 
respectively, and Vay tne shearing strain. dhe Strain components for small 


displacement are 


ou, 
sae Se 
ou 
ou ou 
ee 
Vw YW * & 


In terms of the potential functions, the equations for the stress components 


become 





a 


2 2 a = 
C= ee + 2u (oe. 8 Y) 





oe Jx- ae ae Oxoy 
CO Co oO 
: a9 d¥, dF : 
o =H (2 Soy 52 + : ) (2-7) 
oO Co oO Co 
Lon cms a9 oY 
a ——- 3 + — + “— wa 


In polar coordinates 


VN 23 F Os 


or 


oo, 12,139 _ 123% 
a os = 30° of ae 
(2-8) 
ov Le, ey 1 ey 
are TOF LP ag? on ious 
Displacements in polar coordinates become 
z 1 ov 
at 3 "> 3 | 
(2-9) 
_1@_ 
"=F ee 


where u and v represent radial and tangential components of displacement, 


respectively. 


otress components may be written as 


Co Co Co 
_ ag 1 oY 1 3 a] 139, lw. oy 
or as oe or ye st] r 2 mae 00 «Ora 
= ea a. oy 1 av 1ay (2-10 
eos pe oe 20 2s eon 2-10) 
5° oC oC oO 
- ee pean lov ay 3 loy 1 oy 
"9@ = yf ee 3e at at ee r orod — we ay 
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Wiere Oo ae eae represent respectively the radial, shear and circum- 
ferential (hoop) stresses in the medium. 

In the following development of the equations, the incoming wave, 
considered a step function of time, is either a dilatational or shear wave. 
The step wave was chosen as a convenient form from which, through the principle 
of superposition using Duhamel's integral, the effects of any time dependent 
wave may be approximated. It is shown by Kolsky (4) that for either the 
incoming dilatational or shear wave, the reflected wave from a plane body must 


have both dilatational and shear components in order to satisfy boundary con- 


ditions. The same holds in the present case. 


2.11 Dilatational Wave 

The orientation of the incoming dilatational (compressional) wave 
with respect to the spatial coordinates and time are shown in Fig. 1. Owing 
to symmetry about the x axis, the radial stress, hoop stress, and the radial 
component of displacement are even functions of 9, while the shear stress and 
tangential component of displacement are odd functions of @. 

In the expressions previously given, the dilatational potential may 


now be separated into two parts 
as Pin a Pout 


where 


a se % (acre, t) 


represents the potential of the incoming step wave, and Dee represents the 
potential of the outgoing dilatational wave. 

The displacement and velocity of a particle behind the plane wave 
front can be derived from the potential of the incoming wave, with notation 


Gien dige deal sh io. 





The velocities in 


a a, %' (xre,t) 


polar coordinates become 


cos 9 


Oo 
pea 
Oo 
2 sin 9 
pc, 


13 


(2-11) 


(2-12) 


The radial and tangential displacements at the boundary r = R due to the 


incoming plane stress wave are expanded in Fourier series as a function of 


time as 


Ue ee) 


Oo 


» a_(R,t) cos 


ti 


II 


viR,t,¢) = » b_{R,t) sin 
Neat 


and the velocities as 


OO 


u({R,t,¢) = » a(R, t) cos 


n=O 


CO 


¥(R,t,@) = : b_(R,t) sin 


The Fourier coefficients are determined as a function of 9 


n=1 


ng 


ng 


ng 


ng 


(2-13) 


(2-14) 


» @ measure of the 


degree of envelopment of the shell by the incoming stress wave (Fig. 1), by 


evaluation.et tHe 


integrals 
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0 
a 
CO 
a = cos 8 cos nO dé 
n TPC. 0 
Q 
2G i 
b= ke f Sane sla me wae 
n TOC, 0 


The indicated integrations are performed using the orthogonality properties 


of the sine and cosine functions to get the following result 





Sin 2 al = 16) 
oO sin 25, 
ae ee mats Qo o+-—————— WS oe 
n(in) mpc, i 2 
sin(n-1)9, sin(n+1)0, 
n=1 + 1 a ee 
(eri) 
; sin 2c, 
Oo he 2 ie 
e 4 p a 
aan) ~ «foe 
ea sin(n-1)@, sin(n+1)@, 
n-1 - n+1 n= y 9; 


The subscript (in) has been added to designate the coefficients due to the 
incident wave. Co is a time dependent variable which varies from O to tm 
during the transit of the incoming wave across the cavity, after which it 
becomes a constant equal to @. 

The stresses in the medium behind the wave front are as shown in 
Fig. 1, where ae a ratio relating the stress in a direction perpendicular to 
the direction of wave propagation, is derived from the assumption that there 


is no strain behind the wave front in the direction parallel to the wave front. 


It is related to Poisson's ratio, vy, of the medium in the following way 


* The bar over the symbol does not indicate a vector quantity. 








1) 


Fo 2-16 
Vea (25) 
The stresses behind the wave front are in polar coordinates 
= (eon + Sn aey 
ee Pp 
iy 
= a 20 
Orns ( 3 ein 
SA, = =5 (ee ay, eee) (2-17) 
ee Pp 
Stresses in the medium as a function of radius and time are now 


expressed in terms 


0 (rst,%) 


0 .g(rst,%) 


Sog(r,t,@) 


and the coefficients determined in the same 


velocities are 


0 


| et ek 
eaten Om 
ee 
n(in) On 


of Fourier series 


CO 


» A (r,t) cos 


=o 


co 


2 B (r,t) sin 


na 


(2-18) 


oO 


» C (r,t) cos 


=O 


manner as used previously for 


- ee 
(1+v)@,, + ( 5 cin 20. n = 0 
— ee l-v\ .. 
s og = 
(1 ve, + (l+y) sin oa oe) Seen ‘l n=2 
— = = a) 
(1-v)sin({n-2)0, (1-v}sin(n+2)@, A(l+v) sin no, 
n-2 n+2 n 
n=) eee 
sin 4o, 
0 ee Nea 
sin(n-2)é sin(n+2)é 
ald il 
n-2 - hye Dee 
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a l-v) ., 7 
(1+v)@, = i 5 ) ein 204 n=0 
ar a : aa 

= = = 2 = is i g : Reset ien? ] 0 = 
n(in) a (v ve + (l+v) sin 2 1 + ( i Nees aay 1 n e 

(v~1)sin(n-2)@, (y-1)sin(n+2)@, 2(1+y)sin nd, 

st eS 

n-2 n+2 n 


n= 15 ee 
where the subscript (in) again refers to the incident stress wave. 
As noted earlier, the reflected and diffracted wave is represented 
as a diverging or outgoing cylindrical wave from a line source whose origin is 
the axis of the cylinder, and includes both dilatational and shear components. 


The outgoing potential functions are expressed as infinite series, each term 


of which will henceforth be called a mode 


CO 


i 7 £_(r,t) cos né 


D = 
n=O 

Wee = ) g (r,t) Sin me (2-20) 
nel 


It is shown in the appendix that the coefficients are of the form 


CO 





i = (G1)e [ F’(t - — cosh u jecoch au. d 
n n C mn spare 
C1 O uf 
(2-21) 
_ (10) [ n iS 
ae G(t - — cosh u,) cosh nu, du, 
Co 0 a 


n 

where F (t - — cosh u,) and Gt = cosh u_) represciy the ath order 
Cy ad Cy 2 

derivative of the respective functions, n being equal to the mode number 


considered. 


17 


To express the velocity components of the medium at the boundary 
due to the outgoing waves as in equations eo by Svevwo li fenentia tes 
equations (2-9) with respect to time and substitutes the expressions given 
above for the potential functions. This results in the following for the 


velocity coefficients 





oo 
n 
: (-1) n+e 
orca) SY aca f F (61) cosh vy cosh nu, du, 
Cy 0 
oo 
n 
(ay) f n+2 
ei ; G (By sinh U, sinh nu, du,, 


(2-22) 


oo 
Nn 
eel - rre(e) sinh u, sinh nu, du 
0 


n 
i {ay Glee ) cosh u. cosh nu. du 
- 0 oC 


where (ene and (6,,) represent the arguments of the respective functions. 
The subscript (out) refers to the outgoing wave. 

The stresses in the medium due to the outgoing waves are also 
expressed in series form as in equations (2-18), the coefficients being 


found by substitution into equations (2-10) 


(-1)"p 
A ————— 
n(out) n+2 
= 


n+2 » 2 
[ PF (En) cosh Sea 2 cosh u, | du, 


n+2 


Co 


a oO 
(-1) UL f n+e2 : : 
- G (€ ) SiuieciessSinh. nu. au 
0 2 2 2 2 


(2-23) 





aS 


(ab LL n+2 ( . 
Eco oe pe Sinh eu) sinh nu, du, 


oo 
(-1) [ g't(e ) cosh 2u, cosh nu. du 
7 n+e 2 OD 2 2 
eG O 
CO 
(cae) ve n+2 nN 5 
“a(eue a Bar [> (65 ) cosh Sal i. 2 sinh u,] du, 
alk 
oO 
i (-1) nt (e joe ees inh nud 


2.12 Shear Wave 

The orientation of the incoming plane shear wave with respect to 
the spatial coordinates and time are shown in Fig. 2. In this case the shear 
stress and tangential component of displacement must be even functions of 9 
while the radial stress, hoop stress, and radial component of displacement 
are odd functions of @. 

The incoming wave expressed in terms of a shear potential function 
aL 

ln g K(x + et) (2-2h) 

from which the velocity of a particle behind the wave front in polar 
coordinates may be derived as 


o 
1 2] ee 
2 


G (2-25) 
° S 
V =— cos @ 


P&5 


The stresses in the medium behind the wave front are as shown in Fig. 2. In 


terms of polar coordinates they are 





19 


Oo =o sin 20 
fs S 
= Q - 
0.9 G_ cos 2 (2-26) 
= - j QO 
Fp6 of Sin 2 


Velocities and stresses due to the incoming plane stress wave as functions 
of radial distance and time are expanded in terms of Fourier series ina 
manner similar to that of the dilatational wave, except that because of the 
difference in symmetry, the sine and cosine terms are interchanged. The 
coefficients are also found in a similar manner and are given here. ‘The 


coefficients for the velocity series are 


sin 20 
ae ad 
Oo 1 2 
4 So 
n(in) ~ apc 





x 
2 sin(n-1)0, sin(n+1)0@, 





n-1 - n+] Ds pee 
(2-27) 
sin e, ig, = 19) 
. 0 sin 20. 
Dn (in) 7 HPC, "i : yo oe 
sin(n-1)@, sin(n+1)@, 
n-1 r n+1 N= 2,),0°- 
and the coefficients for the stress series are 
i Q 
a Saris 1 ; 
a. 1 yy ee 
n(in) — ay * 
sin(n-2)@, sin(n+2)@, 
n-2 Sao ee Me =e dl a) Se ey ee 
i 0 
sin 2 1 ; 
oC n=0 
ae sin ues 
= eX a = 
n(in) 4 Sa realy n=2 
sin(n-2)@, sin(n+2)@, (2-28) 
ten n= 1,5,4,5,..- 


n-2 n+2 
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Cn(in) ~ ~ “n(in) 


G1 is a time dependent variable similar to that for the dilatational wave 


except that now it is a function of the position of the incoming shear wave, 


which travels with velocity Cy: 


The outgoing potential functions are expressed modewise as 





cQ 
= i 0 
Pout » f_ (x,t) sin n 
=i 
(2-29) 
Co 
=) 9 
eae » g (r,t) cos n 
He) 
where the coefficients are of the form 
a co 
_ (+1) pe rv 
fete F (+ - — cosh u, ) cosh nu, du, 
ie i 
dl 
(2-30) 
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oo ae 
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The terms have the same significance as in the case of the dilatational wave. 


Modal coefficients are derived as before. The coefficients for the velocity 





terms are 
co 
: rr Gin pute 
= (out ia ee (5, ) cosh u, cosh nu, du, 
C O 
ie 
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eH sinh Us sinh nu, du, 
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and the coefficients for the stress terms are 
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2.2 Equations for the Shelli 


In the following derivations of the equations of motion for the 
shell, it is assumed that the thickness of the shell is small in relation to 
its radius. This assumption permits description of the behavior of the shell 
in terms of its middle surface, and effectively concentrates all of the mass 


in a line thickness. This should be kept in mind when interpreting the results. 
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The effects of shear and rotatory inertia are not considered, since they 
affect only the very high modes, which are not considered in this study. 

The deformation of any point (@) on the shell can then be described 
by specifying the components of displacement (Fig. 3): 


radial displacement, positive outwards 


Uu 
S 


tangential displacement, positive in the direction 
of increasing 9. 


I| 


Vv 
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The strain energy of the shell per unit length can be expressed as 


EI on : BA on ; 
(eo (k)“Ra@ + — e,) Rdg (2-33) 
e 2 0 
0 0 

where Es = "plane strain" modulus of elasticity of the shell 

I = moment of inertia per unit length of shell 

A = area of the cross section per unit length of shell 

k = change in curvature of the cross section 

ea extensional strain in the 9 direction 

R = radius of the cylinder 


Flugge (2) gives the expressions for the curvature and extensional strain 





terms as 
al su, 
eae Boe 
R o@ 
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Ov 


ea aa oe) 
QQ” R ‘9é S 
correct to first order terms. The strain energy equation in terms of the 


components of displacement then becomes 


Ee aI ou, . BA - Ov. 9 
U = 5 ; San + u,) + oR J (<g> + u,) ag (2-3) ) 


ed 


The kinetic energy of the shell can be expressed as 


Om 


ae 
T = , (a, + v2) Rae (2-35) 
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where 


m = mass of shell per unit surface area 


2.21 Dilatational Wave 
In view of the symmetries of the problem, the displacement com- 


ponents of the shell can be expanded in terms of the following Fourier series 


co 


u_(@,t) = » a(t) cos n@ 
n=O 
(2550) 
v (9, +) = » boat) sin n@ 


n=ib 


Equations for the strain and kinetic energies of the shell are then expressed 


as the following quadratic forms in the generalized displacements and 











velocities 
Ein E An Ein — EAn < 
uU = (=—+ ) at eo) oP (an?)? 4 SY +a ye 
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The Lagrangian equations of motion in terms of the generalized displacement 
coefficients a_ and b_ are 
sn sn 


d ,or OU 
ae (3 la 0a n 
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OU (2-38) 
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where @. and Q,, are generalized forces corresponding to the displacement 


terms a__ and bow respectively. The following equations of motion result 
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Now to determine the generalized forces om and Q the principle 


of virtual work is applied. Consider a virtual displacement corresponding 


to an increment oa, of coordinate a 


Virtual work, by definition = Q., : oe 


The external forces acting on the shell are 
oO 
G -) i cos ng 
18 n 
n=0 


oO 
Oo = » B sin ng 
rg n 
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where AL and B represent the sum of the forces on the shell due to the 


incoming and outgoing stress waves 
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A (in) - A_ (out) 


dj 
I 


7 B_(in) + B (out) 


The coefficients ‘oe and B. have been discussed previously. The work due to 


a virtual change of the generalized coordinate oe 
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Virtual Work = A een nO + 8a 
0 n Sia) 
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The generalized force term is determined by equating the two expressions of 


work, from which we get 
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and by a Similar computation 
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2.22 Shear Wave 
Because of the difference in symmetry associated with the shear 
wave as compared to the dilatational wave, displacement components in the 


case of the incoming shear wave are written as 


ui(9,t) = » a Mt) sin ne 
leet 
a (2-41) 
od 
v(9,%) = ). vg q(t) cos ng 
n=O 


Expressions for the strain and kinetic energies in terms of the generalized 


coordinates 7 and os now become 
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The equations of motion determined in a manner similar to that of the 


dilatational wave are then 





EI ee EA A. 
@ on \ (=n ) : ri : Co (2 . nb.) i. mm 
mR 
(2-43) 
nk A B 
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2.24 Effect of Additional Mass 

The equations of motion derived thus far assumed no mass within the 
shell. The simplest way to consider the effect of additional mass is to 
assume that it is distributed within the shell symmetrically about the axis. 
The total quantity can be assumed equal to 2nRm', with m' the magnitude of 
the added mass measured in terms of the surface area of the shell. If it is 
assumed to move with the same velocity as the mass center of the shell, its 


kinetic energy can be expressed as 





ON Co 
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JU a 7: (eos COs oe wes La 6)ae (a-lh) 
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for the case of the incoming dilatational wave. if U, and v, are replaced by 
their Fourier expansions, and if the indicated integration is carried out, 


the kinetic energy T' due to the additional mass may be represented as 


ee ae oe a 
Le wake eee 1 see 





This kinetic energy must be added to the kinetic energy of the shell which 


was derived earlier. Since the above expression affects only the n = 1 mode, 
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the other modes need not be considered. With this additional kinetic energy, 


the equations of motion for the n = 1 mode become 
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If now we add the above equations, we get 
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By making the appropriate substitutions for a 1 and boy the etuations of 
motion can be written as 
oo (a eb) = Sm 
ak Ro en el ~~ mi 2(m+m' m 2(m+m' 
= _ = (2-6) 
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ee Ro = Silat 2(m+m' m 2(m+m' 


Note that these equations are similar in form to the equations derived earlier 
without the additional mass; the effect of the additional mass merely alters 
the right-hand side. 


In the case of the incoming shear wave, the kinetic energy of the 


adqational mass 1s 
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| . WATE : : ‘ 
= ame (a. sin 6 + Vv, cos Q) 20 (Cele) 


and by a similar method, the equations of motion for the n = 1 mode are 


written 
E A B 
i" oi ( SD yt 2th ee m' 
eal Ro ol ei em 2(m+m' m 2(m-+n' 
_ - (2-48) 
: BA b ) 7 OL ime i Fi emt+n' 
sl ec sl m 2(m+m' m  2(m+m' 


The significance of the additional mass on the numerical results is discussed 
in a later chapter. Inclusion of a flexible support for the additional mass 
is also possible by changing the above equations; however, this problem is 


not treated here. 


2.3 Boundary Conditions 

It is assumed that the shell is attached to the medium permitting 
no differential displacements at the boundary between the two. Thus, con- 
tinuity of stresses and displacements are maintained at the boundary. A 
convenient way of satisfying the boundary condition of equal displacements 
of the shell and medium is to equate the corresponding velocity components 
of the shell and medium. This is done modewise for both the dilatational 
and shear incoming waves by merely equating the coefficients of the velocity 


terms as follows 


0. 
il 


ee a_(in) + a_(out) 


(2-49) 
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2.4 Summary of Equations in Non-Dimensionalized Form 


Tt is convenient to express the equations derived earlier in 
non-dimensionalized form by introducing certain dimensionless variables and 


parameters. Let 


ey 
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er 
sn Ss 
==: = T= 2-50 
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be the non-dimensionalized form of the displacement and time variables. 


Parameters are expressed as 


E 
a 
2” E 
wet 
p Pp 
ie (1+y) (1-2y) 
Vv co 
(2-505) 
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Vt = R 
os 
Ta = R 
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The last two parameters given above become unnecessary when considering an 


unstiffened shell since for this case A and I are given by 


iz _ 1 49 7 
A=t and a (2-52) 


Although the equations were derived considering an unstiffened shell of 
uniform thickness, the response of a shell with thin, closely spaced stiffeners 


can also be approximated by specifying A and I separately from t. The 
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equations of interest in non-dimensionalized form are 


ad = Pl (8, ) = Ql_(F,G) 

B = P2 (8,) + 92 (F,G) 

7 (2-53) 
a, =N1 (@,B.) + P3,(0,) + 3, (F,G) 

B = N2 (a, .B,) i PH (8, ) + Qu (F,G) 


where Pl_, P2, P32» and PH are functions whose values are determined 

directly from the position of the incident wave. Ql» Q2 > Q3 > and QO are 
functions of the outgoing shear and dilatational waves; and Ni and Ne, are 
functions of the displacements. These various functions in non-dimensionalized 
form are derived from equations given earlier. Given here are those for the 


case of the incoming dilatational stress wave. 
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a) Other Equations of Interest 


Dilatational Wave 


The hoop stress in the shell is determined from the relation 


eee 
where 
= 1s ae 
Jaana oo Ss 


By substitution of the Fourier expansion of the displacement components the 
hoop stress can be expressed modewise as a fraction of the absolute amplitude 


of the incoming stress wave, as 


0° 
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09 
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The bending stress is determined from the relation 
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which may also be expressed modewise as 


o_, a eee 
ical = d » ie) @_ cos ng (2-56) 
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where d is the distance from the neutral axis of the shell to its extreme 
fiber. 
Stresses within the medium at any radius r may be determined by 


summation of modal stresses as follows 
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where A Bo and C are coefficients previously described as generalized 


stresses of the incoming and outgoing stress waves. 


Shear Wave 


The corresponding equations for the incoming shear wave are, for 
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where all terms are as defined before. 
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CHAPTER III 


METHOD OF SOLUTION 


5. 2 iGene nad 

In this chapter are presented the numerical techniques used in 
solving the problem. Since the methods employed are similar in theory for 
both the incoming dilatational and shear waves, only the solution to the 
incoming dilatational (compressional) wave will be discussed in detail here. 
The equations (2-53) to be solved consist of two pairs of coupled integro- 
differential equations in the generalized coordinates of the shell and the 
potential functions. The integral terms in these equations contain elements 
of the outgoing wave potentials, both dilatational and shear, whose values 
must be determined through application of the boundary equations at each 
instant of time considered. 

In the numerical solution, time is taken to be the independent 
variable, and is expressed in terms of the half transit time of the incoming 
wave across the cavity, T = te,/R. Increments of time are expressed as a 
fraction of the half transit time; if N is defined as the number of time 
steps for the wave to travel one radius, then i/N denotes the elaspsed time 
after the ith step has been taken (Fig. 4). 


The position angle, G1: can be defined for any step in time as 


9, = are cos (1 - i/N) (3-1) 


After total envelopment, i.e., for i/N > 2, 9, becomes equal to 7. 


af 


4.2 Numerical Integration of the Potential Functions 
Lamb (4) notes that if a point source f(t)dz is located at the 


pointe x =) Oyeye= 10,940 ="7), its effect at a distance r from the origin in 
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the xy plane can be represented by the equation 


N fence 


= ) az 


5 F(t E 


Yar J reere ri 
which may be represented graphically as in Fig. 5. To obtain the effect of 
a line source of density F(t) on the z axis, spherically symmetric point 
sources with the same variation in time can be assumed to be situated all 
along the axis, and their combined effect represented as an integral over all 


the point sources 


oO 
: x Cize 
1 [ F(t BN a2 ) az 
= 00 


9 => =a (3-2) 
an =a N re4ze 


[pemcinttsowol Manus Intnl oy bo PlUS intimruy are shown here in order ce 
represent the general case. It will be shown that for a disturbance beginning 
at some definite time, the limits of the integral may be taken as finite. 

Now let 


renee “(1 + oo /r°) 


2 2 
r cosh Uy 


and the above integral can be written as 


00 


= J F(t ao cosh u,) du, (3-3) 


from which the graphical interpretation in Fig. 6 is made. 

Consider for the present the determination of the integral values 
at the boundary r = R. We define t = O to be that time when the first effects 
of the outgoing stress waves (i.e., from the point source at z = 0) reach the 
boundary. As time increases the effects of other point sources progressively 
further away along the z axis will reach the boundary in the xy plane. Thus 


it can be seen that though the limits of infinity are given in the integrals, 
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for finite times the limits can be represented by finite values. The maximum 
value of the variable R cosh Uy in equation (25) that need be considered can 
therefore be represented as a function of time. At t =0 


R cosh Ur m =n 


that is, Vim is zero for time zero. For any time later 


R cosh Urn hee tc, 


in the case of the outgoing dilatational wave (Fig. 6). R cosh Un locates 
the point farthest along the z axis whose wave front has just reached the 


boundary in the xy plane, Cy being the velocity of wave propagation. In 


non-dimensionalized form this can be written 


a 
cosh Us = ae tS ae v 


where the terms are as previously defined. Since the integration is over all 
point sources whose wave front has reached or passed the boundary in the xy 


plane, the integral expression of equation (4-3) can be written with new 


limits 
“Im 
R 
P= [ F(t - = cosh u,) du. 
O ale 
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cosh u, = 14+ 6 (3~)) 


where € now becomes the variable of integration varying in value from 
Oto t. In terms of € the integral expression becomes 


= 


7 W(=1+7=€ jae ; 
® joe (3-5) 


Note that a singularity of the integrand occurs at € = 0, for all values of Tt. 


Paul (6) has presented a numerical technique for evaluating integrals of this 
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general type taking into account the singular point. The method is essentially 
a modified trapezoidal rule in which weighting coefficients are obtained. If 
the range (0,t) of the integral is subdivided into subranges each of equal 


length At, the above integral for any time T = 14AT can be represented as 


(i-1) (m+1)At 
QP = » [ — (3-6) 
m=O mAtT 26 


where f(€) is a polynomial of first degree in € which approximates the 


quotient 


F(-1+t-¢) 
N1+6/2 


as a linear function over each subrange. The results of the integration of 
equation (3-6) are then put in the form 
(i-1) 


QM = SOI Ne Ps | (3-7) 
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where 
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represent a set of weighting functions which can be computed for any interval 


(3-8) 


being considered. For large values of m these functions were written ina 
form more convenient for computational purposes by asymptotic expansion, to 


get the following result 
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For an illustration of the variation of AM and BM functions with m, see 
eevee Oe R. is termed the multiplying factor which for this particular 


integral is equal to 


All of the various integral expressions for the potential functions were 
simplified to the form of equation (3-6) thereby permitting use of the above 
weighting factors for all integrals. However, the expression for the multi- 
pilyines factor R. varleoeior tne dit ferecmreinverrals. 

The integration for the outgoing shear potential y is handled in a 
similar manner, except for a minor modification which results from the fact 
that the shear wave travels outward with velocity Cos 
Consider for illustrative purposes some shear potential represented 


as 


oO 


R 
Y= J G(t - e cosh un) cosh u,, du, 


Since the incident wave strikes the shell at time t = 0, the shear wave front 
must also reach the boundary at this instant in order to satisfy boundary 


CGOnGierons = aus at tC = 0 


R cosh on R 
and for any time later 
R cosh Yon = het blo 


using the same reasoning as for the outgoing dilatational wave. We can express 


the above equation in non-dimensionalized form as 


cosh ron =] + kt 
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where 
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is a constant relating the velocities. 


Again since u, is a variable of integration we can write 


2 


cosh Uy = 1 + K, € 


where € as in the case of the dilatational wave is the new variable of 
integration with limits 0 to tT. In terms of this variable the integral is 


now written as 


a - 6) Vk, (1 + k §) ac 
t 7 tal (3-110) 
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Note that the shear wave starts from the z axis at time 1 =-1/k, which is 
earlier than for the dilatational wave since it travels at a slower speed 
(with fictitious material within the shell) and must reach the boundary at 
time equal to zero. The numerical integration of equation (3-10) is 
accomplished in a similar method to that which was done for the dilatational 
potential. The result can be expressed as 

a 

ve » AK, , Sn en : a eal os awa | (3-11) 

m=0 
where AM and EM are weighting factors identical to those derived previously, 
and Qn for this particular integral is 
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The expression for a is again dependent on the form of each integral considered. 
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For an illustration of the variation of the F and G functions with time, 


see Figs. 11 and le. 


4.3 Solution of the Basic Equations 

The equations of motion for the shell are solved by an iterative 
procedure known as the Newmark Beta Method (5) with which values for the 
potential functions, and the accelerations, velocities, and displacements 
of the shell are determined. 

In general the method consists of using a step-by-step integration 
technique over successive time intervals assuming a specific variation of 
acceleration during each interval. If we assume that at time t = (i-1)At all 
values of potential functions, accelerations, velocities and displacements 
are known, the method becomes that of determining the corresponding values 
at a time tT = iAt. For an assumed linear variation of acceleration over each 


interval the equations of interest are 
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where n denotes the mode and i the time. Note that in the above equations 


the velocities and displacements are expressed in terms of a and B 


e ° 3 
sak Mies: 


acceleration components which generally are unknown. 
Values for the acceleration components at time Tt = iAT are assumed 
and the velocities and displacements determined from the above equations; then 


by use of the boundary equations to determine the values of the potential 
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functions, values of the acceleration components can be computed from the 
equations of motion. The computed values are then compared with the assumed 
values and if they agree or are within an arbitrarily established tolerance, 


the step is completed; if not, the values of a ; and B ; just computed are 
4 


3 
taken as new assumed values and the cycle repeated until the criterion is 
satisfied. A detailed discussion of the procedure follows. 


The equations of continuity at the boundary at time T = iAT are 


given in equations (2-49) and are rewritten here in the form 


Q. 
| 


A 
= Ellery + J Ql (Fo,F),---Fy; Gy»G,, -.-G,) at 


(3-13) 


T- 
II 


1G 
Pema ee f oop i GorG,,++-G,) ale 
O 


where Ply and Pe. denote functions related to the particle velocities in the 
medium at the boundary due to the incoming stress wave; their values may be 
determined directly at any time. The integral expressions denote the con- 
tribution due to the outgoing waves. All values of F and G are known except 

B and G. which are to be determined from the above equations. These equations 


are now written isolating the unknown values of F and G terms 


TAT 
a4 = Poe as f QL (Foo Fyo++ °F, _33 Go rG),+--G,_4) al 
CNG 


At (3-14) 
‘ J a eee oe ae 
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The integrals with limits At to iAt are evaluated using the numerical 
technique described in a previous section. 


The next step is to assume values for a ; and B ‘? the radial and 
oe 


, 
tangential components of acceleration of the shell, from which velocity and 
displacement components are determined using equations (3-12). With these 
values known, equations (3-14) become effectively two equations with the two 
unknowns , E and G., Toi ceeCanmici De Ce termined . 


The equations of motion are given in equations (2-53) and are here 


written in the form 


IAT 
“Cc = : Q eeo ; © © e 
Gg = NIG, 43 By) + P5,(9) + J Q3 For Fpo++Fys GorG,.---G,) at 
(Ca 
aA 
Bas = Nea, 43 Ba) + Ph (9, ) + [ Qu (PPP; Go»G,,---G,) at 


heoOMewmbew tie acCecleration Components aa and BG are now computed. If 
the computed values agree or are within an arbitrarily established limit of 
the assumed values, the step is completed; otherwise the computed values are 
used as the assumed values for the next cycle of iteration. 

The procedure described above requires that all parameters at T = 0 
be known, including the initial values Ee and Ce dnatial values of velcerncy 
and displacement of the shell are specified to be zero. From the short-time 
approximation presented in Section 4.7 initial values of acceleration and G 


Verma om@rie@ tOectilso be Zero, and the values of BO to be the following 


I] 


7 7 -1 n= 0 
Toe ae (3-16) 
G n+1 
2(-1) n= ee 





Tie wetort-Line approximation aillso indicated that the radial acceleration 


1/2 Bye 


components near tT = O varied as T and the tangential components as T : 





43 


thus to improve the accuracy of the calculations, equations (3-12) for the 


first step in time were modified to 


oe ne 
ear Ea ral 
(3-17) 
2 . 2: 
rile ee a Fae Sa TL 


Without this modification, the resulvs forr the firs, Step cane be breucm: 
into acceptable agreement with the short-time approximation. 

The n = O mode is simplified somewhat because it has no tangential 
component of displacement, and contains only the dilatational component of 


the potential functions. 


3.4 Stresses of the Shell and Medium 

Solution of the equations of motion yields modal values of the 
functions F and G, and the modal acceleration, velocity and displacement 
components of the shell. The stresses in the shell are determined directly 
from the displacements using equations (2-55) and (2-56). 

Stresses in the medium at any radius and time are determined using 


equations (2-57). They are here written in the form 


olr,iac) E ; fs.c0,) + | nat 4 face at (3-18) 
ie Sean O O 


where S69) represents the stress due to the incoming wave and the integral 
terms that due to the outgoing waves. All are functions of both radius 
and time. Consider the stresses in the medium for any radius equal to r at 
time t = iAt = oe 

a (Fig. 7) is now determined from the equation 


0. = arercos 2 - i/N) (3-19) 
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and takes on values from 


hon). =O 
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@ = 7 for t 2 (1 + a) 


The integral terms must be recomputed for each radius and time considered. 


As an example consider the integral 


oO 


© = , F(t - ee! u,) du, 
O aL 


shal (hee = = ~ 1, where kK. = R/r, represent the nondimensionalized time 
ie 


after tT = O required for the outgoing dilatational wave front to reach the 


radius r being considered. Then for 1 < - the integral 1s equal to zero, 


and for t > 1, the maximum value of (r cosh u, ) can be written 
= 1 
recosh Usa =) ee Cy 
or 
= a ! 
cosh Us 1 + kK. 
where T!' = oe represents the time after the wave front has reached the 


radius r. We express cosh u, as 


cosh u, = 1 + kG 


where € is the same variable of integration as was used previously. The 


integral is then expressed in terms of € as 


T t 


| B(-1+t'-¢) Vic, ac 
OP = mee a me ee ee 
: Vat 1 + a 


2 


(3-20) 


The numerical method of integration used is identical to that described 
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earlier for the integral at the boundary, with the result that the same 
weighting factors AM and BM. are applicable. Note, however, that the limits 
of the integral are now O to T'; also the multiplying factor R iene Claligs 


pambreular Inveszral Decomes 





By a similar analysis, the integral of the shear potential 


oO 
. i 
Y = f G(t - a, a up) cosh uy du, 


can be put in the form 


i ane = i Vk. Ges kK o)ag 





Y = [ ar a ne (4221 ) 
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Not a+ a 
where 
oe co ae 
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For the potential vy 


Ae il. 
Pete 
C le 


Une welgitims faccrors AM and EM remain the same but the multiplying factor 
Qn for this particular integral becomes 


J KK, (lL +k) 


“a kke 
eS iL 


1 + 5 
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3.5 Time Dependent Stress Wave 

The problem considered thus far has dealt with an incoming stress 
wave with a step distribution in time. The results obtained for the step 
wave can be used through the application of Duhamel's integral to find values 
for the response of a shell and medium to incoming waves with any time 
variation. As an example, a stress wave which decays exponentially with time 


(Fig. 8) according to the following equation is considered 


< 
< i 6 
a(t) =, (1 - 2) « (3-22) 
p po iS 


oe represents the time at which the stress wave decays to zero and k is a 
parameter which is related to the shape of the curve. Stresses at any time 
equal to iAt can be found by the application of Duhamel's integral, here 


Written in the form 


a do (t) 
o(iAt) = a, (idt) + I a g_,(idt-t) dt (3-23) 


one is defined to be that stress resulting from an incoming step wave of 
amplitude ee in equation (3-22). In the numerical analysis, the time 


dependent wave is approximated by a series of rectangular sections as illu- 


strated in Fig. 8. 


3.6 Description of the Computer Program 
The problem was programmed modewise for a high speed digital 
computer (CDC 1604) using Fortran langue se- anCmconsidering enlly stiles tai 
three modes. Input data consist of the following parameters: ~ 
(1) The time intervals at which computations are to be performed, 


expressed as the number of intervals required for the incoming wave to travel 
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one radius (one-half transit time). The degree of accuracy achieved is 
dependent on this parameter, the smaller the interval the more accurate 

the results. However, the machine time required for a given number of transits 
of the incoming wave varies inversely as the square of the interval size ap- 
proximately, thus some sacrifice in accuracy is necessary to reduce the time 
of computations required. 

(2) The total time over which the computations are to be performed, 
expressed as the total number of time intervals to be considered. Generally 
speaking, all values seemed to have reached their asymptotic (static) values 
within ten transit times of the incoming wave across the cavity. 

(3) The ratio of the moduli of elasticity, E/E, where E is the 
"plane strain” modulus for the shell and E the modulus of elasticity of the 
medium. 

(4) The mass ratio, p./; where P., is the méss of the shell andy, 
the mass of the medium per unit volume. 

(5) Poisson's ratio of the medium. 

(6) The ratio of the thickness of the shell to its radius. When 
considering a shell whose area, A, and moment of inertia, I, are not directly 
related to the thickness, A and I must be specified separately. 

(7) The amount of additional mass within the shell expressed as a 
fraction of the mass per unit surface area of the shell. 

(8) The number of radii to which stresses in the medium are 
desired. 

(9) The time intervals at which output data are desired. 

(10) The angular increment at which output data is to be computed. 


Because of symmetry only values between O and 180 degrees need be considered. 
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Addit readin puemauanlLi lies Gin the case or whe exponentially 
decaying stress wave are 7. and k, parameters which define the duration and 
shape of the stress pulse, respectively. 

Cumemencavea consist of the acceleration, velocity, displacement, 
and stress components of the shell for specified angles and times; and the 


stresses in the medium for specified radii, angles and times. 


4.7 Short-Time Approximation 

We a @heck Om the accuracy of the machine solution for short tames, 
the boundary equations and the equations of motion were solved approximately 
by making a series expansion of all pertinent functions in terms of time as 
the independent variable. Although the following discussion is limited to the 
ease of the incoming dilatational wave, the basic principles are the same for 
either of the types of wave considered. 

The expressions for the velocities and stresses in the medium 
around the boundary due to the incident wave all have been written thus far 


in terms of 9 the position angle of the wave. In terms of nondimensionalized 


1? 
tame t, all functions of a, can be written in terms of t using the following 
relations 

cos a, =l-tT 

cin, = Varfi-zr-d?- 2. (3-24) 

A) = Jor F ct = + ~~ ro ts a 


The integral values representing the effects of the outgoing waves 
are also expressed as functions of t. The following example will illustrate 
the technique used to accomplish this. For example, consider the transforma- 


tion of the integral 
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QP = [ F(t - -- cosh u,) du, 
0 a 


As was shown in Section 4.2, cosh Uy in the anvegral varies in value ter an, 
given time t from 1 tol +t. It is convenient for purposes of analyzing the 


integral to represent the variable cosh u, as follows 


cosh u, = 1 + wt (3-25) 


where t is now a fixed value in the integration and w is defined as the 
variable of integration whose value ranges from 0 to 1. The function 
R 
F(t - — cosh u,) = F(-1 + t(1-w)) 
Cy iL 


is expanded in terms of a power series as 


oO 


F(-1 + 1(1-w)) =F + : 7,7 (1-0)? (3-26) 
ie 


where Ys are unknown coefficients of the series. Since 


the integral can now be written as 


du 





— +Z (wt)* t+. . | dw 


p = “Ef afr, +r) +e] fa SE | ae 


from which after performing the indicated integration 
=N2v |F + (27 - =r; ) a < (Cest (3-27) 
oc aren Coe 


The integrals which represent the effects of the outgoing shear 
wave can also be transformed in somewhat similar manner. Consider for example 


the integral 
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Co 
R 
Vy = f G\Ge= == cosh Up) cosh u, du, 
O 2 
Co tc, 
Tambois cacemeoscneu. varies from 1 to i+ Kk.t wnere kK = == and 1 = =——. 
2 e c Cy R 
Therefore we write 
cosh Us, = ts BU (3-28) 


where W is now the variable of integration ranging in value from O to 1 in 
the integral. The function 
G(t - ee cosh u,) = — + t(1-w)) 
Cy 2 kK. 


is expanded in terms of a power series as 


oO 
aL as i 
o(- + x(a-¥)) = 6, + > eyrta-u) 
tk 


where €, are unknown coefficients of the series. From equation (3-28) 


N2Ok t 
el. = = 
fe oadtu 








Wk T 
a fe + % (uk) : os dw 


The integral written in terms of the variable Ww is now 


Nk 1 . 7 k wt 
y=— —G¢ die | B S ae fa wk dw 
is > ja lce + €, 1 ( Ww) + elas + Wk t 


and performing the integration 





y =VoK + (2 G.) + O( e, + ~ kG) + n?(... | (3-29) 


The same basic technique is applied to all the integrals so that now the 
effects of both the incoming and outgoing waves can be written in terms of T. 
Substitution into the continuity equations and the equations of motion yield 


the following for any mode 
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Bezels 
(3-30) 
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ere fl,, f2,, f3,, and fh are coefficients which contain certain elements 
2 Go ee and €. of the potential functions; @l1 and fe are known functions 
* the displacements. 


To solve the above equations, the displacement components of the 


ell are expressed as Frobenius (8) type series 


~ 
ey 
Qa = aU 
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(3-31) 
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i=l 
ere Dj» Qs; and c are unknown coefficients. Substitution into equations 


-30) results in the following set of equations 
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oo oo oO ils 
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The coefficient c is now determined by inspection. Then through a step by 
step process which involves the equating of coefficients of like powers of T, 
values of Ps» Qs; ae and €, are determined. The following equations are then 


used to find values of the potential functions, and the displacement components 


of the shell. 


oO 
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F(t) = fe > y,7 
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G(t) = G+ » E,T 
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(3-33) 
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Velocity and acceleration components may be determined by differentiation of 


the above expressions for the displacements. 
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CHAPTER IV 


DISCUSSION OF RESULTS 


4.1 General 

Results of computations performed to determine the effect of the 
various parameters are discussed in this chapter. 

Although equations presented throughout the study have been written 
to include an infinite number of modes, the greater part of the actual calcu- 
lations performed and presented here are the results obtained considering only 
the modes n = O, 1, and 2. It is important to note that during envelopment of 
the shell by the plane stress wave, a Fourier series representation of the 
incoming wave is objectionable in that the series at this stage is slowly con- 
vergent, thus necessitating a large number of modes to accurately represent the 
plane wave. However, after passage of the wave across the cavity, the Fourier 
expansion of the incoming stresses around the boundary results in coefficients 
of all modes except n = O and 2 becoming identically equal to zero for the 
plane dilatational wave, and coefficients of all modes except n = 2 becoming 
identically equal to zero for the plane shear wave. Therefore, stresses due 
to the outgoing waves in modes corresponding to those of the incoming wave 
whose coefficients become zero must also eventually vanish at long times. The 
limited study conducted for modes greater than n = 2 indicated that the 
maximum effect of the higher modes occurs within one transit time of the 
incident wave and rapidly decays, thus contributing relatively little to the 
maximum response of the shell which occurs after several transit times. 
However, for determining the early time response of the shell and mediun, 
the higher modes are significant and should be considered in further extension 


Gite pairs. Work. 
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In the tables and figures to be discussed, quantities given in 
non-dimensionalized units are defined by equations (2-36) and (2-50). Stresses 
are given in units of the absolute value of the amplitude (lo,| or }o_|) of the 


incident wave; a negative stress means a compressive response to an incoming 





compressional (Fig. 1) or a positive shear wave (Fig. 2). The physical 
pmepernies Cf tic shell relative vo those of the medium, as well as wie naire 
of the incoming wave are indicated on the graphs. Unless otherwise stated, 
the shell is considered to be an unstiffened one so that its cross sectional 
area and moment of inertia are related to the thickness as given by equations 
(2-52). Except where indicated, there is assumed to be no additional mass 
within the shell. Numeral subscripts denote the mode number. 

The shell and medium have been assumed to exhibit linearly elastic 
beiavier throughout thelr Suress hasvories, whch for the practical propillem 
does not permit evaluation of any spalling or non-elastic effects. 

Values of stresses given are in addition to those which exist prior 
to the arrival of the incident wave. For the elastic case, the effect of 
prior stresses such as those resulting from the overburden may be taken into 
account by merely adding them to stresses caused by the incident wave. 

Honmclam ty Ol presentation Sndebecause Of the impractuicabilli gaa 
tciudine solutions for all possible permutations of the parameters invelyeca, 


the discussion in this chapter is limited to a few representative cases. 


4.2 Modal Response of the Shell and Medium 

Figures ll and 12 illustrate the shape of the modal components of 
the dilatational and shear potentials obtained in the solution to a typical 
problem. It appears that a singularity occurs at one transit time in the 


case of the F functions resulting in the slight irregularity of the curves at 
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thie point. olnce computed valucs of stresses were determined to be rather 
insensitive to relatively large variations in values of the F and G functions, 
the effect of the irregularity would seem to be slight. 

Figures 13 through 20 show modal acceleration, velocity, displace-~- 
ment and stress components for the shell and stress components for the medium 
at the boundary, as they vary with time. Static values shown were computed 
using the method given in Appendix B. 

The high accelerations computed near the beginning are not truly 
representative of the actual case, since they are the result of assumptions 
made earlier in deriving the equations for the shell. The shell was repre- 
sented by a line describing its middle surface which permits no variation in 
accelerations, velocities, and displacements of particles through the actual 
thickness. Also, no provision was made for refraction of the incident wave 
tarough the shell dinins.” These Wimitatsons restrict the applicability of 
the solutions to a shell whose thickness is small relative to its radius. 

The n = 1 mode is primarily a translational one which accounts for 
the rigid body translation of the shell after it has been enveloped by the 
incident wave. Thus, it can be seen that the velocity components for this 
mode approach constant values equal to the velocity of the medium behind the 
incident wave front, and displacements grow without bound reaching a straight 
line variation with time. Note that the stresses contributed by this mode 
reach their peak values within one transit time and quickly damp out, 
approaching zero asymptotically. For the incident shear wave, the n = O mode 
is also a rigid body movement which accounts for rigid body rotation, and 
contributes little to the stresses. 

Modal quantities obtained are coefficients of Fourier series; 


therefore, the total response or effect is determined by adding the coefficients 
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Mubeipiltedspy whelappropriate sine and cosine verms for any desired angle. 
Figures 21 through 24 show the time variation of stresses in the shell and 
medium for various angles when the first three modes are summed. The maxi- 
mum stress in the shell for any angle may be determined by adding the bending 
stress to the hoop stress. This is indicated in Figs. 21 and 24 by the 
dotted line above the hoop stress. 

Figure 25 is given to illustrate the relative magnitudes of the 
hoop stresses in the shell and medium for several thicknesses of shell. This 
also shows the effect of varying the relative thickness of the shell on the 
hoop stress in the medium. The dotted line indicating the hoop stress in 
the medium for an unlined cavity was obtained from the report by Paul (6). 

Figures 26 and 27 show how the relative thickness of the shell 
affects the radial and shear stresses in the medium. 

As was discussed earlier, stresses in the medium for any radius 
can be determined by reevaluating the integral terms which represent the 
effects of the outgoing waves, and adding them to the Fourier expansion of 
the incident wave. Figure 28 shows the modal and total radial, hoop, and 
shear stresses which were computed for a time equal to 10 transit times. 
These values are compared later with the static stresses, but on this figure 
the static stresses do not differ by more than the thickness of the lines, 
and therefore are not shown. The time variation of the radial and hoop 
stresses in the medium for various radii are shown in Fig. 29 for the incident 


dilatational wave, and in Fig. 450 for the incident shear wave. 


4.3 Short-Time and Asymptotic Comparisons 


A method for obtaining a solution to the problem which is accurate 


for very short times (t << 1) was presented in Section 3.7. This was desirable 
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to validate the machine solution and to determine the effect of the size of 
time interval selected. The results for a representative problem are shown 

in Figs. 31 through 33 for the incident dilatational wave, and in Figs. 34 
through 36 for the incident shear wave. Four terms of the series representing 
the F and G functions, and three terms for other quantities were used in the 
short-time solution. 

As can be seen from the graphs, good agreement was obtained for 
very small values of time, somewhat shorter time being obtained for the shear 
wave as compared to the dilatational wave. The shorter time results from the 
nature of the forcing functions (Eqs. 2-19 and 2-32) which indicate a more 
rapid rise in the incoming stresses for the incident shear wave. 

Within the range of time for which the short-time solution is valid, 
decreasing the size of time interval for each step of the machine solution 
resulted in closer agreement between the two methods, as is to be expected. 

It also indicated that the stresses and displacements are not as sensitive to 
variations in the interval size as are the F and G functions. 

At the other end of the time scale, i.e., at a relatively long time 
after passage of the incident wave front across the cavity, another check on 
the accuracy of the machine solution is afforded by the asymptotic approach 
of all values to the static results. Figures presented thus far have shown 
that the static condition is approached well within ten transit times. 

The time interval used in the machine solution affects the stability 
of the results for long times. This is indicated in Fig. 37 which shows the 
variations in computed values of the displacement components for mode n = 2 
at relatively long times, for different time intervals. N represents the 
number of time steps required for the incident wave to travel a distance equal 


to the radius of the opening, and 1/N defines the interval size. Note the 
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smaller graph which shows the percentage difference between the computed 
values at ten transit times and the static values. 

As the interval is decreased the machine solution at long time 
approaches the static solution more closely. Below a certain size of time 
interval, there is little difference in the results, which indicates asymp- 
totic convergence to the correct solution. For this particular problem, 

N ~ 30 seems to be “critical” in that for N < 30, wide variations in computed 
values occur. As the mode number increases, the "critical" value of N in- 
creases rapidly, and the requirements of computer storage and calculation 
time become decisive factors which make impractical the study for long times 
of modes much larger than 2. The interval size selected for all problems 
solved (exclusive of the study to determine the effect of the interval size) 
was set equal to 1/40 (N = 40) of the half transit time of the incident wave. 

Hor ete static case, only the modes n = O and 2 yield values other 
than zero. Table 1 compares values of shell stresses and displacement com- 
ponents in these modes obtained from the computer solution to a particular 
problem at a time equal to 10 transit times, with the static solution. Most 
pairs of values differ by less than one percent. Comparable agreement of 
stresses in the medium at various radii are shown in Table 2. Values in 
Tables 1 and 2 were obtained from the solution to a problem whose parameters 


were: =O aa 15 00 Ny = +05, and v = .25. 


Up 0 

Although the figures and tables presented above were for a particular 
problem, the discussion given is applicable to all problems which were solved. 
Changing the physical characteristics of the shell and medium within the range 
of values studied had hardly any effect on the degree to which the long time 
machine solution and the static solution agreed. Also, the absolute value of 


all quantities which should asymptotically approach zero became less than 


-O0005 in each case well before ten transit times. 
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hk Effect of Parameters 
Studies were conducted to determine the effect of each of the 


following parameters: 


_ “8 _"s 
'k ~ E hae 
1, =< vy =v (medium) 
t  R 


for an unstiffened shell without additional mass. Since it was impractical 
to take into account all permutations of the above parameters, a basic shell 


where 


Ip = 4e6 ie BO 


Me = oD V = ite 


was considered from which each parameter was separately varied to determine 
its effect on the resulting stresses and displacements. Calculations 
performed were only for the case of the incident dilatational (P) wave. 

Of particular interest was the determination of the maximum dynamic 
stresses and displacements (not including rigid body translation) due to the 
incident stress wave. Tables 4 through 6 compare the maximum values obtained 
in the machine solution with the static solution. "DLF,"” termed the dynamic 
load factor, is defined as the factor by which the displacement or stress 
produced by or applied as a static load should be multiplied in order to obtain 
the maximum dynamic value. Figures 38 through 41 are given to graphically 
illustrate the variations in stresses and displacements, both static and 
dynamic, in the range of parameters considered. Stresses are given in units 
of io | and displacements in non-dimensionalized units defined by equations 


(2-36) and (2-50). 
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Table 4 and Fig. 38 illustrate the effect of increasing the 
relative thickness of the shell. The range of values selected for ue is 
probably much greater than is practical or justified by the assumptions of 
the analysis, and was considered only to determine the trend of the results. 
As Te Homie reascdadioplacements and Hoop Stresses of both the shell and 
medium decrease; however, the bending stress of the shell and the radial and 
shear stresses in the medium at the boundary increase. Figure 25 shows the 
time variation of the hoop stress in the shell and medium for various thick- 
nesses of shell, including the case of the unlined cavity. 

Increasing Te Results 4anoa rapid) increase in the displacements and 
stresses of the shell and a much lower rate of attenuation of the hoop stress 
in the medium. See Fig. 439 and Table 4. 

Table 5 and Fig. 4O show that the parameter i has no effect on the 
static results but does affect the dynamic response. Increasing ie increased 
the maximum response in both the medium and shell. 

Figure 41 and Table 6 show that within the range of Poisson's ratio 
for the medium considered, as vy increases displacements and stresses in the 
shell decrease, with very little additional reduction of the medium hoop 
stress. 

Additional mass within the shell was shown to contribute relatively 
little to the overall response. This is mainly due to the simplifying assump- 
iGioms waach were made in deriving the equations of motion to account for the 
additional mass. Figure 42 shows the effect of the presence of additional 
mass equal to 407 times the mass of the shell on the dynamic response of a 
particular shell. Additional mass decreases the displacements and increases 


the dynamic stresses in mode n = l. 
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4.5 Response to Time Varying Incident Wave 


The results obtained from a solution to an incident wave with a step 
variation in time was shown to be useful by application of Duhamel's integral in 
determining the response to any time varying stress wave. The case of the 
exponentially decaying wave was considered, the results of which are 
illustrated in Fig. 43. Substantial reductions in the maximum stresses can 
be expected as the duration becomes smaller. In Fig. 44 is shown the effect 
of a linear rise in the amplitude of the incident wave. Note that for a 
wave with a linear rise followed by a step variation in time, very little 


decrease in maximum stresses occurred. 


4.6 Comparison with Previous Work 

Baron (1), using a different method of analysis, investigated the 
dynamic response of two shells with different physical characteristics sub- 
jected to an incident plane dilatational wave. Figures 45 through 49 compare 
his results for the modal values of stresses and displacements of the thin 
shell, with results obtained by the method of solution given in this study. 
Similar results were obtained considering his so-called stiff shell. Although 
the shape of the response curves can be considered similar in the two reports, 
the magnitude of the dynamic response in Ref. (1) seems consistently higher 
than in the present report, and the long time results asymptotically approach 
values higher than the static solution. The following table compares the 
maximum modal stresses, for Poisson's ratio of 0.25,obtained in the two reports, 
and also shows values for the static case. Stresses are given in units of 
lo, The bending stresses for the stiffened shell cannot be compared directly 
Since d, the distance from the neutral axis of the shell to its extreme fiber, 


Momnetestaved in Retr. (1). 
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Thin Shell otiff Shell 











Maximum Stresses 


Ref. Current 
ay Work 


Maximum Stresses 





Ref. Current 
(1) Work Static 





Quant. Mode 


a9 O -4.79 ~4.08 ~4.67 
Gog 2 5.42 5.20 2695 
Oo O 0.55 4 08d 4.67¢@ 
5 ° 1 Seoul 11.894 





4.7 Conclusions 

Conclusions drawn from the results of the analysis are: 

(1) The method which has been presented is practical for effective 
computation of the dynamic response of a cylindrical shell embedded in an 
elastic medium when subjected to plane dilatational or shear waves. The 
solution presented herein is believed correct since it was checked by independent 
methods at short and long times. 

(2) Peak stresses and displacements in both the medium and shell 
occur sometime after the transit of the incident wave across the cavity; of 
the problems solved, both the average and mean time at which the peak values 
occurred was equal to 3 transit times. The dynamic effect measured in terms 
of the ratio of the maximum stresses and displacements to the static values 
varied within a relatively small range. The average value of this ratio was 
ite and tie mean, il.l1l. Thus, for the cases considered at least, the maximum 
stresses and displacements to be expected for a particular situation can be 
roughly approximated by determining the static values and by multiplying them 


iomoetaetor of i.l1. 
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(3) The largest stresses in the shell occurred for a relatively 
thin liner in a medium with a low modulus of elasticity, low mass density, 
and low Poisson's ratio. The greatest reduction of the hoop stress in the 
medium as compared to the unlined cavity results from a relatively thick 
liner in a medium with low modulus of elasticity and high mass density. 
Additional mass within the shell has relatively small effect on the dynamic 
Sugecees. 

(4) The practical value of tunnel linings to reduce the maximum 
stresses in the medium depends on the several conditions mentioned above, and 
on the magnitude of the incident stress wave. Under favorable conditions, 
reduction of stresses on the order of 30 percent or more is possible. However, 
for materials such as granite, smaller reductions can be expected for steel 
liners of practical size. 

(5) Certain assumptions made concerning the behavior of the shell 
have limited the applicability of the analysis to relatively thin liners. 
Future studies of the behavior of thick shells would indicate the effect of 
the approximations used herein. The analysis has also been based on the 
assumption that the behavior of both the liner and medium is linearly elastic 
at all times. Perhaps a more desirable condition would be one in which some 
inelastic behavior is permitted to take place in the medium surrounding the 
shell or in the shell, or one in which some inelastic energy absorbing medium 
such as cinders or foamed plastic surrounds the shell. As a subject of future 
study, it is recommended that the behavior of thick shells, and thin shells 


surrounded by some energy absorbing layer, be considered. 
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APPENDIX A 


DERIVATION OF THE EXPRESSIONS INVOLVING THE POTENTIAL FUNCTIONS 


A.l General Form of the Potential Function 


The dilatational potential function must satisfy the wave 


equation 
de,1%,1 39 _1 3g (4-1) 
oC r or Co 2 2 2 - 
or xr oO Cy ot 


m in the case of the incoming dilatational wave is an even function of 9 and 
thus can be expressed as a cosine Fourier series 


QQ 


= : f_ (r,t) cos ng 


(A-2) 
n=O 


where Hn is a function of r and t representing the modal coefficient of the 


potential function. By substitution of the above expression into equation 


(A-1) we obtain the equation that must be satisfied by f_. 


Nae of 2 ei 

pepe cock tor eee (A-3) 
pee e Dn areas 22 

or r Cy ot 


The general solution of f is assumed to be of the form 


n 
a =r R (r,t) 


(A-k) 
iiien R, 15 some function of r and © satisfying 
a 2 
3 a 2 f° an (2n+1) OR, 
eh Ty or (A-5) 
ot OME 


ee Ro is a solution it can be shown that the corresponding equation for 


Roel is satisfied by 


al OR, 


n+l ror 


(A-6) 
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By repeated application of this result it can be shown that equation Ss) 


ioe cauisi Led Dy 
- (A-7) 
where Ro is the solution of 


S°R S°R OR 
Ont of O A i | (A-8) 
se 1 yr? ror 


Tie soluvion to this equation for the case of a wave diverging from a center 


is 
oO 


R_ = f F(t - — cosh u,) du (A-9) 
0 C1 ale i 


Hor the preet Of this, see Lamb (4). Therefore, the coefficient to of the 


potential function is written as 
n : 1g 
f =r (==) [ F(t - — cosh u_) du (A-10) 
0 Co 1 Al 


Hememodes ©, 1, and 2, this can be written 


oO 


- is mS) du, 


Kh 
lI 


oO 
Pee an (ea eeosht urdu 
L Cy 0 ale IL ile 


Fh 
II 


co co 
ale i 12 2 iL , ' 
D 5 [ EF (€,) cosh Uy du, + a [ EF (oy cosh Uy du, 
Cy O al O 


where primes indicate the derivative of the function F with respect to its 
argument. Note that the expression for the n = 2 mode contains derivatives 
of two different orders. It is convenient to express the function ro in 


terms of derivatives of a single order which may be done through integration 
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by parts. For the n = 2 mode we then wish to change the second integral to 
an integral involving the second derivative of the function F. Integrating 
by parts gives 
1m “Ln <Lm 
[ gdq = En 2 [ nds 
O O O 


where 


E = Hes) 7 sinh Uy 


r tt e = 
deca: = Cc, F (6, ) sinh Uy du, dy = eosh Uy du, 


Thus 


oO 


a . a elm 
i ia , 10 : 1 . 
[ F a) cosh Uy du, = e f F (6,) sinh u, du, + E (6, ) sinh aD 


The second term can be shown to equal zero since for its lower limit 
sinh uy = O; and for its upper limit, F' is the integral of the function F- 
at the wave front which under the assumed initial conditions, does not exist. 


Therefore, the coefficient of the potential function for the n = 2 mode can 


be written as 


co 


_ IL : ae ae : 
i ae [ ee =~ cosh u,) cosh 2u, du, eine 
Cc) O 1 


Pao cin llareprocess tier cOciiiIclents Of The potential function for eny sore. 


mode can be reduced to the form 


co 


n 
if 3 cu f F(t - — cosh u,) cosh nu, du (A-12) 
n n c ie 1 ile 

Cy O iL 


The shear potential function must satisfy the wave equation 


il 
2" ror re Sor 
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in the case of the incoming dilatational wave must be an odd function of @ 


and can be represented as 


00 


y = » g_ (rt) sin né 


iat 


and proceeding exactly as in the case of the dilatational potential the 


general expression for the coefficient a is found to be 





A.2 Velocity Terms 


Come a eee 
g = G (t - — cosh u,) cosh nu, du 
; C, 2 e 


, (A-14) 


The equations for the velocity components of a particle on the 


boundary due to the outgoing waves 


< c 
Il lI 
[Fh OA OY 
S 
She 
OI Mm ct 
ct a 
' Pit 
OY OY 
3 ND Sle 
<— 
ct ct 


The components of velocity and the 


as 


OO 
u = » a (r,t) cos 


n=O 


oO 


y= » (Ge, 2) 


new 


w/ 
- 
=) 


ao 


aa » £ (r,t) cos 


lel 


oO 


Y = » g_(r,t) sin 


Valter 


ane 


(A-15) 


potential functions are expanded in series 


ng 


ng 


(A-16) 


ng 


ng 
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Substitution into equations (A-15)}) gives 


2 of nN 0g. 
eee Sue 
aft aa 
if g 
° n n n 


The expressions for the coefficients fo and g, are given in equations (A-12) 
and {A-14). Using these in the above equations and by the application of 


integration by parts we get the following velocity terms 








co 
a (out) = - come | pote ye ) cosh u, cosh nu, du 
n 7 antl 0 nek 1 ak 
L 
co 
(1 eee ee 
* nel ; G (64) sinh Uy sinh nu, du, 
ze 
, (A-18) 
° : ( =o i [ Ta ee ‘ ; 
b_ (out) =e orl ; i (6) Sinh Uy sinh nuy du, 
ut 
oO 
Eas nae 
rf nel f c ee cosh u, cosh nu, du, 
2 


A.3 Stress Terms 
Otresses in the medium due to the outgoing waves are written in 


series form as 


oO 


oe , A (r,t) cos n@ 


oO — 
n=O 
oO 

re a » B (rst) Sin n@ (A-19) 
n=l 


oO 
» C (r,t) cos n@ 


n=O 


o@ 





(O 


By substitution into equations (2-10: the coefficients As B and C. are 


expressed in terms of the coefficients of the potential functions, ri and a 








af og ng 
A {out} = (A + 2u) 5" eee ee | 
rh rvor 2 
3 r 
5 as : me = no een 36, 
ir Yr r &n ca a or 
en er See (A-20) 
ea ee r or 2 “n 2 5 
or r 
of 2 de - 
c (out) = mas | a ee 2 g = Qo ‘a -n | 
rl r i oe ies yall ee nN or 
sf Og ng 
+ x Et nm Be ut 
5 2 ror 2 
ae ke 


Again, using the equations for oa and go given earlier we can write the above 


coefficients finally as 





OO 
n 
A (out) = (=) : Bona ) cosh nu ES ecen ieee 
n 0 Liu an ib 


ale forte ae 
as | sinh eu, sinh nu, du, 


B bout ) move fe baa 6.) sinh 2u, sinh nu, du 


one al i 
(A~21) 
(ee (aye fa h h 2 
ne cosh nu, cos oe du, 


oO 
a ee r 2 
C _fout) emerge . F (6, ) cosh nu,j~ = 2 sinh'u, | du, 


S a 
L 


& 
eeu ri ante re Binweeus Sinha nu. du 
; 0 2 ieee 
Za 
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APPENDIX B 


STATIC SOLUTION 


B.1 Dilatational Wave 

The static solution presented here is based on the application of 
the theory of stress functions presented by Timoshenko (9). Under static 
conditions, it can be assumed that at large distances from the boundary of 
the cavity the state of stress in the medium is equal to the stress field in 


the medium behind the front. In polar coordinates this is 


a 
= Q 
0 - (——) on - (= v) g,, cos 2 
= = Gs ieee (B-1) 
ro 2 Pp 
oe aya ae va a cos 26 
09 


ue is the stress in the medium in the direction of wave propagation and vo 
is the stress parallel to the wave front. V> derived from the assumption 


ceiac there 1S no strain parallel te the wave front, is equal to 
voz (B-2) 


It can be seen from equations (B-1) that the n = O and 2 modes describe 
exactly the free field stresses. Therefore, the unknown stresses at the edge 


of the cavity can likewise be expressed in terms of these modes 


- P - P_ cos 28 


O 
vet O Co 


(B-3) 


j 9 
oy, 5 Sin 2 


0 
where Po: P, and 5, are the unknown modal components of stress acting at 


Ppiewpeulldary ps uhese Same Stresses must act on Che shell. Thus the boundary 


stresses in the medium and shell can be illustrated as shown in Fig. 9. 





Te 


The equations of equilibrium in the medium under plane strain 


conditions are satisfied by the following expressions for the components of 





stress 
1m, 1 92 
Cnr = mom) 62 ae 
r o@ 
2 _1 99 1939 (B-4) 
r¢é 200 Yr oroé@ 
Peon 
sss 
“For 
where 2 is the stress function in terms of r and @. 
Eel m= 0 Mode 
Lie wenecral solmuaen of the stress function for the n = © modem. 
given by Timoshenko as 
e 
=e lom a Mas 
from which can be derived the stress components 
k 
et a 
a 
0 a (B-5) 
k 
S59 =- at eM 
a 


mie coefficients K 


emo or = KR 


M= 


For the conditions 


and M are determined from the states of stress atr—~ ~ 


Ro a es P| 


- (=) o 


(B-6) 
P 


of plane strain, the strains are 





where 


(9 


SE leg hema ie 
7 E , ate x 06 
A — 
: | agg = 5 oe (B-7) 
+g 
Vie 
E . 
5 modulus of plane strain 
l-v for tne medium 


Displacement components are found by suitable integration of the following 


equations for the 


eo 


strains in the medium 


_ su 
- or 
vu ldo 
= (B-5) 
_~is,wiy 
~ y¥ 90 or r 


Displacements in the medium for the zero mode are 


a 
R 


Re 
R 
The corresponding 


compressive force 


Eda 


where 


Pe] (B-9) 


displacement components of the shell with an exterior 


Po are found in Flugge (2) to be 


POR 

eee 
A 

Be (B-10) 
=) 

ae 
= T. ua modulus of plane strain 

S for the shell 
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The unknown stress P is now determined by equating the radial displacement 


components, to get 


@ee eG 
P = (B-11) 
5 = 
ae a le 
ie . 
Ss 


After P. is determined, stresses in the medium at any radius can-be found 
using equations (B-5). Equation (B-10) gives the displacement of the shell; 


and the hoop and bending stresses of the shell are, respectively 





Eee 
Cag = - 
A 
(B-12) 
0 
sB__ ust 
oy eR 


B.12 n= 2 Mode 


The general soiution in terms of a stress function 2 is given as 


a 


Q = x + ae "| cos 20 
° r 


from which the stress components in the medium become 


6M oN 
ra = = ax on + = cos 20 
- ie ic 
6M 
og = E: E a cos 20 (B-13) 
EM = ON . 
C9 = ax ~ vi - a sin 20 


Coefficients are determined as before from the states of stress at r-—-~ wo 


aad Cr = R 


(B-1}) 





2 


fee 


Reel : _ 

5 a P, 28, | 
ee — 
a | --¥) oD ais P, ats s, 


of 


A 
I 


The displacement components of the medium at the boundary can be expressed 








as 
— 2(2-y)s 
U 1 a (5-v) J 
See een eae ee P. + ———— || cos 20 
R E a B) e 3 
(ee) 
eel 2(2-v) (5-v) | 
a =n lata ) Be 5 S| sin 20 
and the corresponding components for the shell are 
Ys Por? %o | 
es | = = — - = —| COs 20 
R |. Eg YL 
Ss S 
(B-16) 


2 3 S 2 
2 R R 2 pak R ‘ 
ee + eRe AD] ein 


a 
a 
R 


The unknown boundary stresses P. and 5 are determined by equating the dis- 


2 
placement of the shell and medium at the boundary. The hoop and bending 


stresses in the shell are then found from the following equations 





PR 25,R 
G99 = (Gr + aR) Cos 29 
B-17) 
tR? ( 
ari Ey “- S,) cos 20 


Displacements of the shell and stresses in the medium can be determined using 
equations given earlier in the discussion. Total static stresses and dis- 


placements are merely the sum of the modal values for any angle. 
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B.2 Shear Wave 


The free field shear wave stresses expressed in polar coordinates 


include only the n = 2 mode 
OG =O Meo 20 
eng S 
= Q - 
gg = 9, cos 2 (B-18) 
S99 = -9. Sin 20 


where om is the amplitude of the incoming wave. The static solution is 
obtained exactly as in the case of the dilatational wave except for the 
interchange of sines and cosines resulting from the difference in geometry. 
The resulting equations are 


a. Displacement components 





u i R? S R? 
Ss o 2 . 
or ——— | sin 20 
cee 1 Bor 
; (B-19) 
a (2, Ry es ee 
Ro. jig ~SL A Tm 121 A 
oo 3B. 
be. Shell Stresses 
PLR 25_R 
Tag = - es ah | sin 20 
3 (B-20) 
eee (2P. + S_) sin 20 
sB 121 DOs 
ec. Medium Stresses 
CCD en 
OR (2A + =e a) sin 20 
0 r 
0 ~~ (on -  _ 2d) , 26 (B-21) 
rd Be Ye oe 2 


CC noes 
a 





io 


where 


os 
as 


of 


Pik vaca, 
J Fo = 8 


In the machine program, the n = O mode is purely rotational and 


a 
an 


thus we get a static value of tangential displacement. This is calculated 


as follows 


© 
II 
Oi 
= 
{| 
PO's 
Q 


— 
(B-22) 


I 
eb 
I 
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Q 
he 





78 


DILATATIONAL WAVE 


a Dev. From Dev. From 
Mode Quant. Stat. Mach. Static Mode Quant Stat. Mach. Static 


0 -aQ 4.000 4.006 15% -O% 4.613 4.599 30% 
0 -Oep 4.000 4.006 5 B 3.845 3.856 229 


0 Osb 100 .100 00 Oep 3-077 3.124 1.20 





Osb 346 345 .29 


* Stresses given in units of IGpl 
Displacements defined by Eqs. (2-50) and (2-36) 


SEEARSWANE yime = 20°R/co 


Mode Quant. Stat. Mach. 
(3 5.000 4.996 
ou 13.840 13.930 
3 11.535 11.539 


-Op, 9.232 9.149 





Osb 1.038 1.045 


TABLE 1 


COMPARISON OF STATIC AND MACHINE VALUES FOR MODAL RESPONSE OF THE SHELL 
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DILATATIONAL WAVE time = 20 R/c) 





2.0 R 3.0 R 4.0 R 

Stat. Mach. | Stat. Mach. 
615 .616 638 .639 
.718 .720 696 .697 
226 .229| .271 «273 
384 .386 -363 366 
340) =. 342 





* Stresses given in units of |G! 


SHEAR WAVE | time = 20 R/c2 


1.0 R 2.0 R 3.0 R 4.0 R 
Mode | Quant. | Stat. Mach. | Stat. Mach. | Stat. Mach. | Stat. Mach. 
2 Ovy | -.460 -.4551 1332 -331 | .675 .678| .814 .815 
2 Ove 921 .913 | 1.283 1.286} 1.151 1.154 | 1.090 1.092 





2 ~ Jeg 2.616 2.636 | 1.101 1.104 | 1.020 1.022 | 1.006 1.009 


TABLE 2 


COMPARISON OF STATIC AND MACHINE VALUES FOR MODAL STRESSES IN THE MEDIUM 








Ove, (90°) 580m (90° 
TeV gtatic Max. pur* | tal || static Max. DIF | ta 
Ol -9.21 -10.12 1.10 3.22 -2.49 -2.73 1.09 3.13 
05 -7.08 -7.7% 1.09 2.88  -2.01 -2.19 1.09 2.88 
1005.56 == 6.15. 1.102.250 2067-08 1. 5 
m5) e600 = See) 03:25 0CiCwDWHS SC=65 S'S 
-20 0-394 = 52K 2252302650015 2 25 
ee el | Tire ma (90" 
| Static Sta 
01 lll 3.36 -. 











* Ratio of Max. to Static # Mme in transit times at which Max. occurs 
TABLE 3 


EFFECT OF RELATIVE THICKNESS OF SHELL 
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Odes —, e Seew(90°) 
= ic ae 








2 -4.16 -4.57  #+.1.10 2.76 -2.27 -2.51 1.11 2.76 
4 - 7.08 -7.74% 1.09 2.88 -2.01 -2.19 1.09 2.88 
6 + 9.33 -10.11 1.08 2.76 -1.81 -1.97 1.08 2.76 
10 -12.57 -13.41 1.07 2.76 -1.55 -1.65 1.07 2.76 
16 -15.72 -16.47 1.05 2.88 ~1.31 -1.37 1.05 2.76 
[awa] ww | w [oe YR] oe | 

‘ we [sate [mes | oe | 
2 23 26 1.12 3.25 - .21 ~ .23 1.13 3.12 
4 45 50 1.11 3.25 - .35 - 39 +11. 2.88 
6 64 mil peell 93212 - 47 - .51 1.10 2.76 
10 1.01 1.12- 1.10 3.12 - .63 - 68 1.08 2.62 
16 1.53 1.68 1.09 See eco - .83 1.06 2.62 





Static Max. 
pe 


- 4.60 - 5.12 
- 8.61 - 9.0 
-12.22  -13.37 
-18.65  +-20.20 
-27.11 -29.15 





"Ratio of Max. to Static FT me in transit times at which Max. occurs 
TABLE 4& 


EFFECT OF RATIO OF MODULI OF ELASTICITY 
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Ooo, Q0 Cee. an° 


"p Static Static 


-7.08 -7.67 
-7.74 
-7.97 
-8.30 
-8.70 





fae [S| we | S| oe Le 


1 -8.61 -9.38 1.09 3.85 

3 i -9.530 1.10 3.00 : 4.339 «1.13 3.50 
5 . -9.64 1.12 2.88 . 44d LoS 3.62 
T i -9.77 1.1% 3.00 : 4.45 1.16 2.76 
10 i “9.95 1.16 2.62 : 4.72 1.23 3.00 


*Ratio of Max. to Static fTime in transit times at which Max. occurs 
TABLE 5 


EFFECT OF MASS RATIO 
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Ooo, 907 O80m ( xy 
fstacse | ms | nut | tof | statte | axe | our 
- 8.33 - 9.17 
ol - 7.96 - 8.69 


sé ai! T42 = 8.08 
025 = 7.08 = 7.7% 
a4 - 5.76 - 6.75 


Osb (90° 
wa 








Orr ~" 00° 
Static Max. tm 





*Ratio of Max. to Static Foime in transit times at which Max. occurs 


TABLE 6 


EFFECT OF POISSON'S RATIO 





8h 


Oy Gre / 


"ay 
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Plane Dilatational Wave Front at time t =0 
Stresses Behind Wave Front Shown in Both Rectangular 
and Polar Coordinates 





Plane Dilatational Wave Front at time t70 


FIGURE 1 
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Plane Shear Wave Front at time t = 0 
Stresses Behind Wave Front Shown in Both Rectangular 
and Polar Coordinates 





Plane Shear Wave Front at time t>0 


FIGURE 2 








FIG. 3 DISPLACEMENT COMPONENTS OF THE SHELL 


Wave front = 
— 





FIG. 4 DEFINITION OF POSITION ANGLE, = 
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xy plane 


FIG. 5 GEOMETRY OF A POINT SOURCE ON THE z AXIS 
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FIG. 6 GEOMETRY OF MANY POINT SOURCES 
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FIG. 7 POSITION ANGLE, ©, , FOR ANY RADIUS 
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FIG. 8 TIME DEPENDENT STRESS WAVE 
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FIG. 9 STATIC STRESSES ACTING ON THE MEDIUM AND SHELL 
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HALF TRANSIT TIMES, Y= tc)/R 
FIG. 11 FUNCTIONS F(-1*7) and G(-1/k,* T ) 
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FIG. 33 COMPARISON OF SHORT-TIME RESULTS FOR INCIDENT P WAVE 
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FIG. 34 COMPARISON OF SHORT-TIME RESULTS FOR INCIDENT S WAVE 
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FIG. 38 EFFECT OF RELATIVE THICKNESS OF SHELL 
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